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We consider the excitonic insulator state (often associated with electronic ferroelectricity) , which 
arises on the phase diagram of an extended spinless Falicov-Kimball model (FKM) at half-filling. 
Within the Hartree-Fock approach, we calculate the spectrum of low-energy collective excitations 
in this state up to second order in the narrow-band hopping and/or hybridisation. This allows to 
probe the mean-field stability of the excitonic insulator. The latter is found to be unstable when the 
case of the pure FKM (no hybridisation with a fully localised band) is approached. The instability is 
due to the presence of another, lower-lying ground state and not to the degeneracy of the excitonic 
phase in the pure FKM. The excitonic phase, however, may be stabilised further away from the 
pure FKM limit. In this case, the low-energy excitation spectrum contains new information about 
the properties of the excitonic condensate (including the strongly suppressed critical temperature). 



PACS numbers: 71.10.Fd, 71.28.-fd, 71.35.-y, 71.10.Hf 
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I. INTRODUCTION 

The Falicov-Kimball model (FKM)^ was introduced 
more than forty years ago with an objective to describe 
metal-insulator transitions and mixed-valence phenom- 
ena in certain hexaborides and oxides. It attracted 
an impressive amount of work within both solid state 
theory^ and mathematical physics-* , dealing with all as- 
pects of the ground state properties and the phase dia- 
gram. Early on, various extensions of the original model 
were proposed^'^ in order to more adequately describe 
the physics and/or phenomenology of specific systems. In 
recent years, much attention has been paid to the pos- 
sibility of electronic ferroelectricity in an appropriately 
extended FaUcov-Kimball model (EFKM)^. This phe- 
nomenon is intimately related to the notion of excitonic 
insulator, which was introduced earlier in a somewhat 
different contexls^r— (see also Refs. ITl| - [l3l ). Finding an 
actual electronic ferroelectric might be important tech- 
nologically, and the experimental search is ongoingi^. In 
addition, it has been suggested that the EFKM might be 
relevant for some systems of current interest, including 
manganatesi^ and URu2Si2^. 

Various approximate methods^^^™ (as well as numeri- 
cal techniques2ii^) were used to determine the phase di- 
agrams of both pure and extended FKM and the carrier 
spectral properties in the ground state. Yet it appears 
that the present article reports the first systematic study 
of the low-energy collective excitations in the excitonic 
insulator (or excitonic condensate) state of the EFKM. 
Even though restricted by the Hartree-Fock mean field 
approach, the results shed light on the nature of much 
discussed instability of the excitonic state in the limit of 
pure FKM. In the case of EFKM, we obtain the exci- 
tation spectra in the excitonic state, and thereby iden- 
tify the underlying energy scales. Importantly, this also 



allows to probe the mean-field stability of the uniform 
excitonic insulator. Until now, this latter issue has been 
addressed only variationally, by considering a finite set 
of competing mean field ground states. 

The spinless Falicov-Kimball model proper involves 
fermions di and Ci in the localised and itinerant bands, 
interacting via a Coulomb repulsion U on-site: 
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where is the bare energy of the localised band. We 
choose the units where the hopping amplitude t and the 
period a of the (d-dimensional hypercubic) lattice are 
equal to unity. Presently, we consider the half-filled (n ~ 
1) case at zero temperature. 

Historically, the raison d'etre of the Falicov-Kimball 
model is the presence of an intermediate valence regime, 
whereby the two original bands are partially filled. It 
is well-known^'^ that within a broad range of parameter 
values in this regime, the lowest-energy spatially uniform 
mean-field solution is characterised by an interaction- 
induced spontaneous hybridisation, A — {c\di). Hartree- 
Fock decoupling of the equations of motion for the 
Green's functions {{cj:ci)) and ({d^ct)) yields 
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for the hybridisation and 
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for the localised-band occupancy rid- Here, N is the 
number of the lattice sites and = E^d — eg with 
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the renormalised relative energy of the locaUsed band, 
Erd = Ed + U{\ — 2nd), and the tight-bindmg dispersion, 

This mean-field solution has its filled and empty 
quasiparticle bands separated by a hybridisation gap of 
2?7|A|, corresponding to an excitonic insulator. If the 
two original bands have opposite parity and ReA ^ 
0, it has a spontaneous dipole moment (electronic 
ferroelectricity)^'^'^. However, the uniform excitonic in- 
sulator state in the pure FKM is unstabl^'^'*, which has 
largely been attributed-:^ to fluctuations due to the local 
continuous degeneracy. The latter, in turn, is associated 
with the phase of each of the operators di in Eq. ([T]) , and 
hence with the phase of (c|di) on-site. 

This strong degeneracy can be broken, and the stabil- 
ity of the excitonic insulator may eventually be restored, 
by any of the terms in the perturbation, 

5n= d\d, -I- ^0 ^ cjd. - -| ^ (c\d, + c]d) - 

(ij) i (ij) 

- Y E { - • ^} (^'^J- - 4^^) + H-C- , (4) 

where t' is the d-band hopping and Vb, bare on-site hy- 
bridisation. Vi (V2) is the spatially-even (odd) nearest- 
neighbour hybridisation, as appropriate for the case 
where the two original bands have the same (oppo- 
site) parity. Ri is the radius- vector of a site i, and 

aS = X]q=i ^a, sum of Cartesian unit vectors. 

When the Falicov-Kimball model is extended by 
Eq.Q, it ceases to be exactly soluble. Identifying its 



ground states and constructing the phase diagram consti- 
tute the subject of a broad on-going research effor t^^'^^ . 
In the present article, we address the issue of stability of 
the excitonic insulating state with a uniform real A, and 
its low-energy properties at the mean-field level (ignoring 
higher-order correlations). 



II. LOW-ENERGY COLLECTIVE 
EXCITATIONS IN THE EXCITONIC PHASE 

We consider the particle-hole excitations of the exci- 
tonic insulator state, which in general are given by 

k 

+Fc{K q) cicu? + 9) 4rffe+ J , (5) 

where the four quantities Fi are arbitrary functions of k 
and q. The energies uj and the wave functions of collective 
modes satisfy the secular equation, 

[x^,n + sn] = LjqX^, (6) 

where the label "eff" implies Hartree-Fock decoupling 
on the 1. h. s.. Upon substituting Eq. ([5]) into we 
find a system of four equations for the functions Fi. For 
example, collecting the terms with ctdg_|__. and those with 

dtc^_i_ - yields respectively 



(w - 4, - t'ej:^^)F+{k, q) + (UA* ~ V~,^.)F,ik, q) - {UA* - y^)F,(fc, q) ^ U [A^ + A^iq)] , (7) 

(w + + t'ej:)F. {k, q) - {UA - V^)F,{k, q) + {U A - V^*^^)F4k, q) = U [-Aa{q) + A(g)] . (8) 

Similarly, 

{UA - V^*^^)F+ {k, q) - {UA* - V^)F^ {k, q) + {uj + Cfe^- - ij^)F,{k, q) = UAM , (9) 

-{UA-V:^)F+{k,q) + {UA* ^V~^^.)F^{k,q) + {u-t'ej:^- + t'e^)Fd{k,q^ = UAM- (10) 
Here, the quantities Ai originate from the interaction term in Eq.([T]) and obey the self-consistency conditions. 
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is the Fourier component of the bare hybridisation, 

^ ^ a=l 

(depending on the relative parity of the orbitals). The 
tilde accents in Eqs. (fTTI - [T4| reflect the fact that in the 
presence oiSH, the r.h.s. of Eqs. are trivially mod- 

ified, and the corrected expressions [see Appendix, Eqs. 
(|A1HA2|) ] should be used here. Ideally, one should now 
solve Eqs. (IT HTU)) and substitute the explicit expressions 
for Fi{q) into Eqs. (|imi4p . Zeroes of the determinant 
V^ujjq) of the resultant system of four linear homoge- 
neous equations for Ai{q) would then yield the spectrum 
of collective excitations for a given momentum q (cf. Ref . 
1^. Since we are interested in the low-energy excitations 
only, it is possible to follow a simpler route as sketched 
below. 

For the pure FKM with SH = and real A, these 
homogeneous equations for Ai (q) have an cjg- = solution 
with Aa ^ —A and At ~ Ac ~ Ad ^ 0, corresponding 
to22 
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[which the reader can verify by substituting into Eq. ^] . 
The presence of an entire branch of zero-energy excita- 
tions is an expected consequence of the local continuous 
degeneracy as discussed in the Introduction above. For 
the EFKM with sufficiently weak perturbation, Eq. Q, 
this branch acquires a small but finite energy that can 
be found by expanding the equations for Fi and Ai in 
powers of w, Vi, and t' . 



III. LEADING-ORDER INSTABILITY OF THE 
EXCITONIC INSULATOR 

Assuming that A (but not necessarily Ag, see Ref. 0) 
is real, we find to leading order: 

2?(^, q) = DUq) ■ {-^y + Mll{q) ■ D^{q). (17) 

Here, the first term conies from expanding 'D{uj, q) of the 
pure FKM in powers of lo: 

D^{q) = -[Y^{q)f-ZMyM-WM-ZM]\ 
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Since the two quantities 
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vanish in the long- wavelength limit, q — > 0, we find 
Di^{q = 0) < 0. Likewise, one can show that D^{q) 
at q — Qo = {tt, 7r(, tt)} is given by2^ 
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and numerically we find that Di^ remains negative 
throughout the Brillouin zone (BZ). 

To leading order in perturbation, the determinant V 
for the EFKM [see Eq. dTT])] at a; = factorises into a 
product of the 3x3 diagonal minor (corresponding to 
Ab, Ac, and Ad), 

DM - [4 - 2Z2m {[YM? + 4[ro(g)]'} + 

+2 {1 - [ZM?] Yi{q) - SZMYM , (19) 
and the diagonal matrix element corresponding to Aa, 

■V^+ViEd 



t'{d + t^) ^ I 
UA^Nd ^ fe 1 
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even, 



When Vi = 0, the latter quantity (and hence^iSE ijj^^ van- 
ishes at g = 0, corresponding to the U{1) remaining de- 
generacy for i' 7^ that was discussed earlier—. When 
t' — Vq.i = 0, a different Goldstone mode is present at 
q — Qq. It is due to the degeneracy of the excitonic insu- 
lator with respect to assigning arbitrary opposite phases, 
ztip, to (c|di) in a checker-board order while keeping Ud 
and |(c|(ij)| = A unchanged^. This can be verified di- 
rectly by analysing the mean field equations at V2 7^ for 
the 2-sublattice case. In all other cases, the appropriate 
choice of signs (t', VqA, ViEdA < 0) ensures that Mu is 
negative throughout the BZ. 

At (f = 0, the quantity Dq is proportional to the deriva- 
tive, dA/dVo, calculated in the unperturbed case, and 
is therefore negative (see Appendix |B]). The equation 
I?(a;, q) = then yields a positive and hence a stable 
spectrum in the long-wavelength limitj^S. It is, however, 
easy to show that at q~ Qq, 
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(20) 



and is positiv^^^. This gives imaginary uj at q = Qq, 
implying that the excitonic insulator is always unstable 
when the pure FKM limit is approached. 
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Contrary to earlier suggestion o^'^^ , this type of insta- 
bility is due to the presence of a lower lying ground 
state, and not to the local degeneracy of the excitonic 
phase in the FKM. This agrees with those variational 
results^^'^* that suggest that a lower-energy ordered state 
is favoured in the pure FKM. The excitonic insulator is 
then stabilised by a small but finite perturbation via a 
second-order phase transition^^'^^, whereas a more ex- 
otic behaviour— would have been expected should the 
vanishing (cjdi) in the pure FKM case2^ be due to the 
local degeneracy. The latter, on the contrary, masks the 
instability in the unperturbed case (by causing the spec- 
trum to vanish identically). Hence in order to see this 
instability, one has to allow for an infinitesimal pertur- 
bation, as outlined above. 



IV. STABILISATION OF THE EXCITONIC 
PHASE BEYOND THE LEADING ORDER 

When effects of subleading order in perturbation are 
taken into account, the structure of Eq. p7|) is preserved. 
In the following, we will keep the most significant next- 
order (in t' , Vo,i, and V2) corrections: (i) Corrections to 
the determinant Do, which is responsible for the leading- 
order instability. We will include these by substituting in 
Eq. (Ull) Do ^ Do + Di, where Dq is given by Eq. (jl9p 
with corrected expressions and ni used in the defi- 
nitions of Yn and Z„. Discussion of the quantity Di{q) 
is relegated to Appendix ICl (ii) In all terms in Eq. ([TTl) . 
we use corrected values of A and n'^ [obtained when A^: 
and fii are used in Eqs. ([5H21)]- This is to account for 

the sometimes appreciable^ change of A and n'^ caused 
by a weak perturbation. 

Numerical evaluation shows that the quantity Di at 
q = Qo typically is negative, implying that for perturba- 
tions larger than a certain critical value t'^^. or Vi^cr the 
spectrum may become real and the stability of the uni- 
form excitonic phase with respect to the excitations of 
the type of Eq. ([5]) may be restored, as shown in Fig. [1] 
(see also Appendix [P]) . Applying our perturbative treat- 
ment to a finite SH is justifiable whenever the resultant 
critical values of t' or Vi are small. The latter includes 
requirements, |Vb,i| < UA and \Vi\, \t'\ < 1. 

The critical values of Vi and t' as functions of Ed are 
plotted in Fig. [21 Wherever a direct comparison with 
earlier variational and numerical results is possible (the 
t' casei^i^) , our approach yields an increased stabil- 
ity region for the excitonic phase, with the disagreement 
becoming more pronounced as \Ed\ decreases. This may 
be due to the above extrapolation scheme being inexact, 
or to a possible importance of more complex excitations 
[as opposed to Eq. ([5])], or to the presence of a first-order 
transitioni^!^ at small Ed- 



V. DISCUSSION 

We shall now attempt to put our results into a broader 
prospective. Whenever the excitonic insulator behaviour 
is suggested, a distinction must be made as to whether 
the underlying physics is that of FKM [with one of the 
bands being (nearly) flat] or of a generic semimetal (semi- 
conductor) with nested Fermi surfaces2i"-i2 . Our results 
show that the difference is not restricted to the par- 
ent bandstructure, but is also manifest in the observ- 
able properties of the system. Indeed, the presence of an 
entire branch of low-energy excitations is characteristic 
only of the EFKM-based mechanism (cf. Ref. [!)• 

A physical difference between the two realisations of 
excitonic insulator is indeed expected on very general 
groundai^ii^. In the case of a broad-band semimetal 
(semiconductor), the excitonic (induced hybridisation) 
gap is typically much smaller than the bandwidth of ei- 
ther band. The particle-hole pairing (exciton formation) 
is then essentially a Fermi-surface effect, involving only 
electrons and holes from the vicinities of the two Fermi 
surfaces, therefore nesting is required^. In this case, the 
physics of excitonic insulator is largely the same as that of 
a conventional charge-density wave. On the other hand, 
in the case of EFKM (sufficiently close to the pure FKM 
limit), the width of the narrow band is smaller than the 
gap, hence the entire narrow band is in principle acces- 
sible for excitonic pairing and the Fermi surface shape 
is no longer crucial. The situation is then closer to the 
actual intermediate valence pictureiii^^ (especially where 
slower dynamics is concerned). We note that the band- 
structure of most compounds suggested as candidates for 
excitonic insulating state^-''^i3'36^ includes a very narrow 
band corresponding to a slightly perturbed FKM. 

With increasing temperature, the long-range order of 
the excitonic insulator state in a semimetal is lost ow- 
ing to the thermal excitations of electron-hole contin- 
uum. The critical temperature (which has the well- 
known snail-like profile as a function of band overlap^ 
or of the interaction strengtb^^ U) is found from the gap 
equation^. It appears that the current vie w-'^^'^^i^^ is that 
a similar scenario should hold for the EFKM. 

Our results suggest that the behaviour of an 
FKM-based excitonic insulator involves a second-order 
transition'^^ mediated by the low-energy excitations re- 
ported above^. While even for the pure FKM case 
the gap equation, Eq. has a non-zero solu- 

tion at low temperatures, T < Ta~UA{T = 0)/kB, 
the actual critical temperature Tc for a weakly per- 
turbed FKM does not simply follow the single-particle 
gap. Instead, the relevant energy scale is determined 
by the collective excitations, z. e., by the dominant 
term in Eq. resulting in Tc <^ Ta. It is 

easy to see that Tc will be roughly proportional to 

V\t'm-\t'cr\), \V2\^vi-Vicr, V I ^0 |( | ^0 I " V^O,c.) , 

or Vl^iKl^il-l^i.-D- 

Finally, we summarise our conclusions that appear rel- 
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FIG. 1: (colour online) Effect of different perturbations 
on the excitation spectrum of the 2D FKM, Eq.(IT| with 
U — 2, Ed = 0.4. Solid, dashed, and dotted lines cor- 
respond, respectively, to the effects of t' (bottom to top: 
t' = -0.04, -0.0565, -0.07), Vo {Vo = -0.06, -0.074, -0.09), 
and V2 (V2 = 0.115, 0.249, 0.35). As the corresponding values 
of A vary between 0.214 and 0.359, the gap in the single- 
particle spectrum is in all cases larger than uj^ for any q in 
the BZ. 
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FIG. 2: (colour online) Critical values of t' (solid line, see Ap- 
pendixlDlfor further details). Vo (dashed), Vi (dashed-dotted) 
and V2 (dotted) in a 2D EFKM with U = 0.5 (weak coupling 
regime, right panel), (7=1 (intermediate coupling, second 
right), (7 = 2 (second left), and 17 = 4 (strong coupling, left 
panel). Results for (7 = 1 suggest reentrant behaviour in Vi,2. 



evant in the context of experimental search for an exci- 
tonic insulator with a narrow (several meV) and a broad 

{W - 0.1 eV) bare carrier bands, and U/W - 0.2. (i) 
in this case would not exceed several tens K, and could 
be much smaller, (ii) A peculiar branch of excitations 
("excitonic phonon"^, either acoustic or optical, largely 
undamped throughout the BZ), will be present in the in- 
frared or even microwave range at T < Tc. (iii) Since 
these excitations involve mostly the narrow-band charge 
density waves, Eq. p^ . one can also expect to find static 
or dynamic charge and/or orbital inhomogeneities in a 
region above T^. While these predictions are based on a 
Hartree-Fock study of a spinless model, qualitatively they 
are expected to persist in a more realistic treatment. 
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Appendix A: Corrected mean-field equations 

In the presence of perturbation, Eq. (4), Eqs. (2-3) 
hold to leading order only. The exact mean field equa- 
tions read 



k ' k 



for the hybridisation and 



k 



(Al) 



-Yhf, nf=i- 



2 2^iCj: + t'e-,r+A\UA~V.*\^ 

(A2) 

for the occupancy of the (almost) localised band. 

Appendix B: Sign of Do at g ^ 

Suppose that A (7^ 0) and Ud solve the mean field 
equations, Eqs. (HIE]), for the unperturbed FKM, Eq. 
([1]). Adding a weak on-site hybridisation Vo [i.e., a per- 
turbation Sn, given by Eq. dH) with t' = Vi:2 = 0] 
results in a small change of the mean- field parameters: 
A — >■ A -|- (5A, nd ^ rid + Snd- Keeping the linear-order 
(in SA, Sud, and Vq) terms in Eqs. (|AHIA2[) yields a 
system of two equations. 
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U J N ^ (eH + 4C/2A2)3/2 
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Since in the limit g — (we write instead of {0, 0(, 0)} 
for brevity) Eq. yields 
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-E 

ON ^ 
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2N ^ (^2 + 4(72^2)3/2 

fe k 



U^A^ 



N ^ (^2 + 4[/2A2)3/2 
fe k 



(B3) 
(B4) 
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Eqs. (jBmB2[) can be re-written as 



Appendix C: Expressions for Di 



Vo 



U 



2Yi{0)-6A + AYo{0)-Snd - [1 + 2^(0)] 
This yields 

i6[ro(o)]2 + [i + 2ri(o)]2 



u 



dA 

dVo 



6A _ 

~ {16[ro(0)]2 + 4[yi(0)]2 + 2Yi{0)}U 



(B5) 

This derivative must be negative, as a thermodynamic 
stabihty condition [with the quantity — Vq playing the 
role of an external field, cf. Eq. (j4])]. On the other hand, 
it has the same sign as the quantity DQ{q} at q — ?> 0, 
which (owing to Z12 — > in this limit) is given by 

Doiq^ 0) = 16[ro(0)]' + 4[yi(0)]2 + 2^1 (0) . 



To leading order in perturbation, the quantity Di{q), 
introduced in the main text, has the form 



Vq + VlErd 

UA 



Ho{q) + ViHi{q), even. 



odd, 
(CI) 

(again depending on the relative parity of the two bands) , 
where 



Ho{q) 



Here, 



-(1 - 2na) [a{Y^ + AY^) - 2Y^Z^ - AY^] + 2{YI + 4yo')(8Yi - SaFo + - ^4) + 

+A{YiZi + 2^0)^3 - 2YiZ| + 2(1 + IQY^ + '52Y;^ - 2aYoYi + 4Yi)Z2 + 2(2 - aY^ + 14^1)^2 + 



+2a(12yo2 _ Y^)Zi - 32(Yi - l)yo^i + 2(0^0 - AaY^Y^ + AY^) + |(1 + 2Yi - YiZ2)cri 
+4(^2 + AY^)[(jj, - 2aa2 + (3 + a^)ai\ + A[a + ^2 - aa^ - ^{l - 2nd)]{2Yo + Y^Z{) 
A (^(2nd - 1) - di - 4ri) {AY^ + Yl) - {2AY^ + IW^ + Yi)Z2 - 



(C2) 



-(leYo - 3211,^1 + 10ZiYi)Zi + 1 + 321^2 + + 2(1 - 2nd) 
1 



YiZ^ + 2^0 
A 



-(1 - 2nd) {~YxZ2 - Zl - AaYiZi + (3a2 _^ 8)(y^2 _^ ^y^^ ^ ^Yi - SaYo] - 2(^1^ + AY^) x 
X (Z3 + 32yo) - 16(^1^ - AY^)Zi - SYo{l + 2^1)^2 + 2Zl - IQY^Zj ~ 2(1 + 8Yi)Zi - 
-16ro(l + 2^1) - 8(riZi + 2ro)ai + 6(r2 + AYS){2a<j^ - ^2) - [Y^ + 4ro2)g_|^ . 
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and the quantities 

ZM - 



2 — And Sq!^ eff—2d , . „ 



- 2d 2a 



^ - 8 + Y^iq) - 16Yi{q) - Aa{l + a^)^ + A{3a^ + l)cri - 12aa2 + Aas 



(C6) 
(C7) 



vanish in the long- wavelength limit, q ^ 0, for the unperturbed case. On the r. h. s. of Eqs. (|C2HC4I) . we omitted 
the argument q of the functions Yn{q) and Zn{q)', note that these are computed using the leading-order expressions 
Aj: and nf, Eqs. (EHS]). 

Since the general expression for H2{q) is rather lengthy, here we restrict ourselves to the case oi q = Qq: 



H2{Qo)= -^[{A + a^)h+l2]\ l-2aYo{Qo) 
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FIG. 3: Values of g^r in the t' ^ case for a 2D EFKM with 
(7 = 1 (squares) and U — 0.5 (triangles). The corresponding 
values of t^r a-re plotted in Fig. 2 (solid lines in the two right 
panels). 



It is easy to seess that /i,2 > 0, hence iJ2(Qo) is negative. 



Appendix D: Further details on Figs. 1 and 2 



As the strength of the perturbation, Eq. (4), is de- 
creased towards the corresponding critical value, t'^^. or 
Vi^cr, the collective excitation spectrum softens at a cer- 
tain q — Qcn viz., Lo^^^ = 0, signalling an instability of 
the excitonic insulator. The spectra shown in Fig. 1 
correspond to qcr ~ Qo = {""jTrj. While this is indeed 
a typical case, other situations are possible, especially at 
smaller U. This is illustrated in Fig. [S] In the momentum 
space, the line {0, 0} — {tt, 0} — {tt, tt} — {0, 0} was probed, 
and Qcr was found to lie on the {0, 0} — {tt, 0} — {tt, tt} 
segments. 
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